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Abstract 

We present the results of an age determination study for pre-main sequence stars in the 
Ophiuchus molecular cloud. The ages of eight pre-main sequence stars were estimated from 
surface gravities derived from high-resolution spectroscopy. The average age of the target 
stars was 0.7 Myr. By comparing the individual age and the near-infrared color excess, we 
found that color excess decreases gradually with a constant rate and the lifetime of the inner 
disk was determined to be 1.2 Myr. The estimated lifetime is nearly a half of the time compared 
to that of the pre-main sequence stars in the Taurus molecular cloud estimated with the same 
method. This result indicates that the disk evolution timescale depends on the environment of 
the star-forming region. 
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* Based in part on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 


1 




1 Introduction 


The circumstellar disk of a pre-main sequence (PMS) star is created through the contraction of the 
rotating molecular gas core, where planets form. Disks dissipate due to gas and dust accretion, growth 
of grains, and gas dispersal driven by ultraviolet radiation. Since the formation and the evolution 
of planets are strongly dominated by the disk evolution timescale, an accurate estimation of this 
timescale is desired. 

There have been many investigations to estimate the disk evolution timescale. The early es¬ 
timation was conducted by Haisch, Lada, & Lada (2001) who conducted the JHKL photometry of 
six clusters. They found that 50% of the PMS star disks dissipate in <3 Myr. This relationship was 
refined by Mamajek (2009). The relation between the cluster age and the disk fraction was approx¬ 
imated by an exponential and the e-folding time was estimated as 2.5 Myr. On the other hand. Bell 
et al. (2013) determined the age of the clusters using color-magnitude diagrams, and indicated a disk 
lifetime of 10 — 12 Myr. Many other studies estimated the disk evolution timescale based on pho¬ 
tometric observations (e.g., Hernandez et al. 2008; Muzerolle et al. 2010; Ribas et al. 2014). From 
the low frequency of transitional disk objects, whose disk has an inner hole due to disk evolution, 
compared to PMS stars with optically thick disks, the evolution timescale of the inner disk was esti¬ 
mated to be of the order of 10^ yr (Skrutskie et al. 1990). The disk evolution timescale varies with the 
method used for age determination. Moreover, age determination of a PMS star based on photometric 
observations suffers from uncertainties in distance, extinction, and excesses arising from the accreting 
heated disk. 

The surface gravity, g, estimation of PMS stars is an alternative method to derive the age. The 
surface gravity of the photosphere of PMS stars decrease due to photospheric contraction over time. 
In our previous studies, we established surface gravity diagnostics by calculating the equivalent width 
ratios (EWRs) of Fe (8204.9 A) and Na (8183.3 A, 8194.8 A) absorption lines. EWRs are suitable 
tools to determine the surface gravity of a PMS star, as they are free from the uncertainties arising from 
distance, extinction, and the continuum veiling. On the other hand, we have to take into account the 
possibility that the EWRs of a PMS star with a high accretion rate may contain uncertainties arising 
from line emission filling (e.g., Petrov et al. 2011). Previously, we observed field dwarfs and giants 
of which the surface gravities and effective temperatures (Tes) are well known from photometric 
studies and hipparcos parallaxes (Takagi et al. 2010). The relationship between the EWRs and surface 
gravity was established for stars with spectral types of late-K to early-M. Next, we determined the 
surface gravities of 10 PMS stars in the Taurus molecular cloud using the EWRs. The ages were then 
calculated using the evolutionary model of PMS stars (Siess et al. 2000). By comparing the estimated 
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ages to the near infrared (NIR) eolor exeesses, we found that the NIR eolor excess dissipates within 
several Myr (Takagi et al. 2014, hereafter Paper I). The disk dissipation timescale was estimated from 
a single star forming region. To discuss the diversity of the disk dissipation timescale, we present 
the results of the age determination of PMS stars in the Ophiuchus molecular cloud. Similar to the 
Taurus molecular cloud, this is one of the closest star-forming regions, making it a suitable location 
to estimate the age with high-resolution spectroscopy. 

2 Observations and data reductions 

We conducted high-resolution spectroscopic observations of eight PMS stars on 2013 June 25—27 
with the High Dispersion Spectrograph (HDS; Noguchi et al. 2002) with two 2K x 4K CCD (13.5 
/rm square pixel), mounted on the Subaru Telescope. We selected late-K type PMS stars as targets 
to employ the EWR—log g relation used in the previous studies (Takagi et al. 2010, Paper I). The 
target stars are listed in table 1. All of the targets belong to the p Ophiuchus (p Oph) cloud core 
except V1211 Oph, which is located in the Ophiuchus North molecular cloud. Because a spectrum 
with high signal-to-noise ratio (S/N) is needed (S/N ~ 100) to determine the surface gravity and the 
age, single PMS stars with /-band magnitude of < 14 mag was selected. Since the number of late-iC 
type single PMS stars in the large survey conducted by McClure et al. (2010) is 19, our sample cover 
nearly one third of the known late K-type single PMS stars in this region, and almost all of the bright 
PMS stars suitable for the EWR method are included in present study. The optical elements of the 
HDS were set to an optimum value to obtain the Ee (8204.9 A) and Na (8183.3 A, 8194.8 A) lines 
with resolution of ~60000 (0.024 A pixel“A- The integration time for each object was set to 600 — 
9000 s to achieve an S/N of 100 per pixel. The object which requires exposure time more than 1800 s 
was split into several exposures to avoid cosmic ray contamination. AO-type stars were also observed 
as standard stars for telluric absorption correction. Thorium-argon (Th-Ar) lamp frames were taken 
for wavelength calibration. 

The Image Reduction and Analysis Eacility (IRAE) software package was used for the data 
reduction*. Bias subtraction, flat fielding, cosmic ray rejection, and scattered light subtraction were 
conducted before the spectra were extracted. After wavelength calibration, telluric absorption correc¬ 
tion, and continuum normalization, we measured the equivalent widths of the absorption lines (figure 
1) using the Voigt function. The errors in equivalent widths were calculated from the uncertainties in 
the continuum level. The measured equivalent widths of the target stars are listed in table 1. 


‘ IRAF is distributed by the Nationai Optical Astronomy Observatory. 
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Fig. 1. Spectra of the targets stars. The telluric lines were removed by using the standard star. 


Table 1 . Measured equivalent widths of target stars. 


Star 


ew(A) 



Na8183.3A 

Na 8194.8A 

Fe 8204.9A 

ISO-Oph 6 

017 

n ^ 79+0 007 
u.j/z,_o 013 

»j.u/^_0 004 

ROXsV 

U.DZU_o 028 

n 89 1 ~l"0.033 
U.521 _o 031 

u.iuy_o 008 

DoAr 24 

U.DU ^_0 029 

n 707+0 025 
U./6/_o 022 

n nQ9+0.007 

U.U82_o 007 

DoAr 25 

U.DUJ _0 028 

^•/^D_0.036 

n no^+O.OOO 

U.U8:5_0 008 

DoAr 28 

u.^yD_o 013 

/^_0.024 

0.0791°;°°° 

WSB40 

0.4921°;°=° 

u.Doy _0 016 

0.080l°;°J° 

DoAr 32 

0.6821°;°=° 

n 070+0.031 

U.5/U_o 033 

u.uyo_o 015 

V11210ph 

0.494l°;°i° 

u.o/z_o 016 

U.U89 0.008 
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3 Results 


The surface gravity determinations are conducted by employing the EWR—log g relationship de¬ 
scribed in Paper 1. The relationship between the surface gravity and the Fe (8204.9 A) / Na (8183.3 
A) and Fe (8204.9 A) / Na (8194.8 A) were estimated by observing 21 field giants and four dwarfs 
(Takagi et al. 2010). All of these field stars have of ~4200 K, comparable to late-K type PMS 
stars. The surface gravity of each star was calculated as follows: 


M T L 

log^ = log — + 41og — - log — + log^o, 

IVIq 1q Ijq 


( 1 ) 


where M, T, and L are the mass, effective temperature, and luminosity, respectively. We used the 
parallax in van Feeuwen (2007) and the (/-band magnitude from Tycho-2 catalog (H0g et al. 2000) for 
the luminosity calculations. The Tycho V magnitude was translated to Johnson V magnitude before 
the surface gravity calculation (FSA 1997). The mass of each giant star was determined by comparing 
the location in the HR diagram and the theoretical evolutionary tracks (Fejeune & Schaerer 2001). 
The masses of dwarfs were estimated by using the spectral type - mass relation described in Drilling 
& Fandolt (2000). As mentioned in Paper I, the FWR—log g relationship is expressed as: 
log EWR — logo 


iogg = 


( 2 ) 


where g is measured in units of cm s“k For Fe (8204.9 A) / Na (8183.3 A), the values of a and b 
are 1.522±0.040 and -0.289±0.004, respectively. For Fe (8204.9 A) / Na (8194.8 A), the values are 
1.141±0.024 and -0.294±0.004, respectively. All of the coefficients and their errors are estimated by 
employing the weighted Deming regression^. The error of the surface gravity was calculated with the 
following expression: 


Slogg = 


\ V afelnlO 


Sa ) -f 


'logo — logEWR 
62 


66 + 


{ -^- 

\bEW RliilO 


6EWR 


(3) 


Derived surface gravities are listed in table 2, which range from 3.2 to 3.6 dex. These gravities are 
estimated from FWRs, which show the state of the photosphere and are free from uncertainties in 
distance, extinction and continuum veiling. In addition, the weighted average of log g is 3.42, and the 
weighted average of log g of the four stars with log g <3A2 is 3.29±0.06, and that of the rest four 
star is 3.58±0.06. These points indicate the dispersion of log g and the age spread in a single star 
forming region. 

The age of each target was estimated by comparing the derived surface gravity to the PMS 
evolutionary model described by Siess et al. (2000). For the calculation, the Teff of each target was 
determined by comparing the spectral types quoted from Wilking et al. (2005), Andrews & Williams 


■ These coefficients are updated from Paper I due to the magnitude correction of Tycho-2 catalog and adopting the Deming method. 
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Table 2. Estimated properties of the target stars. 


Star 

Type* 

SpT 

Teff 

logs 

^StELT 

Age 

log -Lstar/ Lq 




(K) 


(AT©) 

(Myr) 


ISO-Oph6 

C 

K4.5 

4470 

o.:)u_o 14 

1 54+0.15 
^•‘^^-0.07 

0 43+0 07 

0 94+0-1'^’ 

ROXs7 

w 

K5.5 

4280 


1-Gy_0.04 

o 

o 

1 + 

o p 

b b 

00 00 

0 79+0.14 
^■^^-0.11 

DoAr24 

c 

K5 

4350 


1 1 9+003 
^•^^-0.00 

0.9lt°;f, 

0 47+0-13 

DoAr25 

c 

K5 

4350 

3 34+0.12 

1 1 7+0.02 
^•^■7_o.oi 

0-82lo.24 

^•^^-0.12 

DoAr28 

TD 

K5 

4350 


1 1 1 -1-0.03 
^•^^-0.00 

i.uj_Q 32 

0 42+0+4 
^•^^-0.08 

WSB40 

C 

K5.5 

4280 

Q aC+O.ll 

■^•■^‘^-0.15 

1 04+005 

l-G^_0.02 



DoAr32 

C 

K6 

4205 

o c'i-|-0.17 
■^•'^^-0.17 

0 90+002 
^•^^-0.00 

0.82l«;« 

0 37+0-19 
-0.16 

V11210ph 

C 

K5 

4350 

T 71 +0-11 

1-^G_q Q5 

0.47l“;« 

0 77+0-10 
^-0.13 


* C, W, and TD correspond to classical T Tauri star (CTTS), weak-lined T Tauri star (WTTS) and transitional disk 
objects quoted from Wilking et al. (2005), Andrews & Williams (2007), and McClure et al. (2010). DoAr 28 was 
classified as TD based on the SED in McClure et al. (2010). 


(2007), and McClure et al. (2010) to the spectral type - relation described by Kenyon & Hartmann 
(1995). The results are listed in table 2. The mass of a target star is estimated by interpolating the 
evolutionary tracks. The values listed in table 2 do not include the errors in the Teg of the targets. 
The estimated ages may include an additional error of ~20%, by assuming a spectral type uncertainty 
of ±1 subclass. This is relatively small compared to the error arising from the S/N of the spectrum 
and the fitting error of EWR—log g. The mass of a PMS star will vary depending on the applied 
evolutionary model (Hillenbrand & White 2004). However, the relation between the age and the 
surface gravity is comparable in PMS stars around 1 Mq. Age estimations based on surface gravity 
is less dependent of the evolutionary models of PMS stars. Figure 2 shows the H-R diagram with the 
plots of the target stars. 

The weighted average age of the targets is 0.7 Myr. Due to the brightness limit of observations, 
the target stars selected in the present study are concentrated at the edge of clouds where the extinction 
is small. The median age of the PMS stars located in this surface region of the p Oph core is estimated 
as 2.1 Myr (Wilking et al. 2005). The age estimated in this work is less than a half of this result. It is 
more comparable to the age estimations of the young stars in the core part of the p Oph cloud, which 
is estimated as 0.3 Myr (e.g., Greene & Meyer 1995). 

To discuss the disk dissipation timescale in the Ophiuchus cloud, we estimated the NIR color 
excess (J — K) arising from the heated disk. The NIR excess of a PMS star arises from the inner 
part of the disk where it is optically thick. A large NIR excess indicates that the hot disk lies near the 
central star. As the inner edge of the disk spreads outward, the NIR excess will decrease due to the lack 
of hot disk. The NIR J and K magnitudes of the targets were taken from the results of the 2MASS 
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Fig. 2. HR diagram of the observed stars. The luminosity of each target were estimated from the comparison of the iog g and the evoiutionary model of PMS 
stars (Siess et al. 2000). The solid lines are evolutionary tracks (Siess et al. 2000) for masses of 0.4 to 1.7 Mq . The red dot-dashed lines and the blue dotted 
lines indicate isochrones and isogravity lines. 


survey (Skrutskie et al. 2006). The intrinsie eolor arising from the photosphere was subtracted by 
quoting the value from Kenyon & Hartmann (1995). The reddening caused by interstellar material 
was removed by adopting the reddening raw (Cohen et al. 1981, Meyer et al. 1997) using the value 
of the extinction {Ay) estimated in previous studies. Extinctions were preferentially quoted from the 
study which conducted extinction corrected spectral energy distribution (SED) fittings mainly based 
on near-infrared observations (e.g., McClure et al. 2010), in order to compare with the result of Paper 
I which also used Ay calculated from SED estimations. Ay was calculated using the extinction law 
of Mathis (1990) when the extinction value in C-band was not presented. The adopted NIR color 
excesses are listed in table 3. 

We derived the disk dissipation timescale in the Ophiuchus molecular cloud by comparing the 
estimated ages with the NIR color excesses (figure 3). The J — K excess decays at a constant rate. 
ROXs 7, which is the only WTTS included among our targets, shows a small excess despite of its 
youth. The presence of this object may indicate that the PMS stars in this regions have different disk 
lifetime, or a rapid disk dissipation occurred in a specific PMS star. On the other hand, atomic lines 
of Na, Ee and V in 8100—8200 A show an asymmetric profile, which imply the presence of cool 
spots or a companion star. This profile was not found in spectra of other targets. The estimated age of 
ROXs 7 may include an additional error from this reason. 

We estimated the typical disk dissipation timescale of this region without ROXs 7. The disk 
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Table 3. The NIR color excesses of the 


target stars. 


Star 

J-K 

(mag) 


Ref.* 

ISO-Oph 6 

0.55±0.03 

2.2 

1 

ROXs 7 

0.10±0.04 

6.0 

2 

DoAr 24 

0.26±0.04 

3.5 

1 

DoAr 25 

0.23±0.03 

3.4 

1 

DoAr 28 

0.15±0.03 

2.3 

1 

WSB40 

0.49±0.03 

4.3 

1 

DoAr 32 

0.18±0.03 

5.9 

3 

V11210ph 

0.43±0.05 

1.1 

4 


* Ay ai'e quoted from; (1) McClure et al. 2010; (2) 
Wilking et al. 2005; (3) Cieza et al. 2010; (4) 
Sartori, Lepine, & Dias 2003. 



Age (Myr) 

Fig. 3. The age - J - K excess relationship of the observed stars. Filled plots are the same as figure 2. Open plots of red circles, green square, and blue 
triangles represent Class II objects, transitional disk objects, and Class III objects of Taurus molecular clouds, respectively, with recalculations of the results 
in Paper I (see the text). The solid line indicates the regression line derived from the plots without ROXs 7. The dashed lines are drawn based on the 1 
rr uncertainty of the regression line al J - K excess = 0. The gray dotted-dashed line shows the relationship for PMS stars in the Taurus molecular cloud 
calculated with the stars with J - K > o. 
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dissipation timescale estimated from the J — K color excess is 1.2±0.1 Myr. Errors were estimated 
from the standard deviations of the weighted regression line coefficients. Because the targets of the 
present work is limited to PMS stars with small extinctions, the disk dissipation timescale of the stars 
in the core part of the cloud could be different. Since PMS stars in the Ophiuchus clouds are young 
(~1 Myr), there is a lack of stars with J — K excess ~ 0. Therefore there is a possibility that the tail 
part of the decline of age-color excess relationship will extend to older ages. Observations of faint 
disk-less PMS stars with next generation telescopes or a near-infrared observations of the PMS stars 
deeply embedded in the molecular cloud can reveal more accurate disk dissipation timescale and also 
the specificity of ROXs 7. 

It should be noted that the age and NIR color excess estimation depends on the employed 
relationships between spectral type, Tes, and intrinsic colors. We also calculated the age and NIR 
color by using the relation of Pecaut & Mamajek (2013) and compared with the result shown in table 
2. The age of each PMS stars becomes 17% older on average, although they are still within the 
range of error. In addition, the J — K excess of each star will be 0.1 mag smaller. Nevertheless, 
the estimated disk dissipation lifetime using this relation is 1.2±0.1 Myr, which is comparable to the 
result mentioned above. 

4 Discussions 

The near-infrared excesses indicate that the disk exists in the inner region of the disk. The decrease of 
the NIR color excess implies the outspreading inner hole. Ribas et al. (2014) estimated the disk evo¬ 
lution timescale by comparing the age of the 22 young clusters and the fraction of the PMS stars with 
near-infrared and/or mid-infrared excesses. The timescales estimated from the shorter wavelength 
and the mid-infrared were 2—3 Myr and 4—6Myr, respectively. They suggest that the dust evolution 
in the inner disk is more rapid than that in the outer part. The estimated timescale in the present study 
indicates more rapid disk clearing in the Ophiuchus molecular cloud. 

The dissipation timescale of PMS star disks in the Taurus molecular cloud estimated with the 
EWR method and the NIR magnitudes from 2MASS was 2.3±0.2 Myr^ (figure 3). The estimated 
timescale in the Ophiuchus molecular cloud in the present study is nearly a half of time compared 
to the timescale in Taurus. The difference in disk dissipation timescale between the Ophiuchus and 
Taurus molecular clouds indicates that the disk evolution timescale will varies between star-forming 

^ This timescale is recalculated by using the observational results of Paper I with the revised EWR-iog g relationship described in section 3. The ratio of 
Sc and Na in the near-infrared iT-band was aiso used for the age estimations of PMS stars in the Taurus cloud. The values of a and b in equation (2) 
for Sc (22057.8A) / Na (22062.4 A) are 2.464±0.05f and -0.289±0.003, respectively. For Sc (22057.8A) / Na (22089.7 A), the values of a and b are 
2.987±0.065 and -0.276±0.004, respectively. 
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regions. 

The gradual disk dissipation with a eonstant rate is eonsistent with disk dispersal triggered 
by the disk wind driven by magnetorotational instability (MRI) -driven turbulenee (Suzuki, Muto, 
& Inutsuka 2010). The vertieal magnetie fields in the disk are amplified by the MRI and the gas 
in the surfaee layer of the disk will flow outward. The eonstant deeline of the NIR exeess ean be 
explained with this wind by assuming that dust grains are well eoupled with the gas. The strength 
of this disk wind is determined by the strength of the magnetie flux. The masses of PMS stars in 
the Taurus and Ophiuehus moleeular elouds used for the disk dissipation timeseale estimations in 
Paper I and the present study are ~1 Mq. Thus, the magnetie flux density of their parent eores is 
eonsidered to be nearly equal. The differenee between the disk dissipation timeseales may imply 
that the fraetions of magnetie field diffusion during star and disk formation are different between the 
Taurus and Ophiuehus moleeular elouds. 

Another possible reason for the differenee in disk dissipation timeseales is the presenee of 
high-mass stars. While the Taurus moleeular eloud is in a statie environment, the Ophiuehus moleeu¬ 
lar eloud is plaeed elose to the Seo-Cen OB assoeiation. The radiation pressure and stellar wind from 
massive stars inerease the density and the magnetie held of the moleeular eloud (Kusune et al. 2015). 
This may suggest that the moleeular eloud eores in Ophiuehus eollapsed with a stronger magnetie flux 
eompared to those in Taurus. Wilking et al. (2005) suggested that the surfaee layer of the Ophiuehus 
moleeular eloud was stripped away by the ultraviolet radiation of massive stars. This radiation may 
have the potential to dissipate the disks, as in the eireumstellar disks in Orion (e.g., O’Dell & Wen 
1994). 

The differenee in metallieity may also be responsible for the unique disk dissipation timeseale 
for eaeh star-forming region. Both observational (e.g., Yasui et al. 2010) and theoretieal studies (e.g., 
Ereolano & Clarke 2010) investigated the ehanges in evolution timeseale eaused by the differenees 
in metallieity, and eoneluded that low metallieity ean result in rapid disk dissipation. On the other 
hand, PMS stars in the Taurus moleeular eloud and the Ophiuehus moleeular eloud have metallieity 
nearly equivalent to the Sun (Padgett 1996). As this study was performed using a hypothetieal value 
of the surfaee gravity, further investigations of the metallieities based on aeeurate surfaee gravity 
and veiling estimations may reveal the eause of the timeseale differenee between the Taurus and 
Ophiuehus moleeular elouds. 

In addition, one may eonsider other faetors for determining the NIR exeess. When the dust 
sublimation radius is smaller than the magnetospherie radius, the inner edge of the disk may be 
eontrolled by the position of eo-rotation and henee the rotation period of the star. On the other hand, 
if the sublimation radius lies beyond the magnetosphere, the position of the inner edge ean be managed 
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by the radiation field. Moreover, the shape and the inclination of the disk can have an impact on the 
NIR color (e.g., Thi, Woitke, & Kamp 2011). Revealing these parameters will contribute to the better 
understanding of the disk evolution and its timescale. 

This research was based on data collected at the Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. This work was supported by JSPS KAKENHI Grant 
Numbers 26103708, 26870507, and 25870124. 
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